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Abstract 



We use ISAJET to perform a detailed comparison of the supersymmetry 
reach of the current Tevatron (100 pb~^) with that of the Main Injector 
(2 fb~^) and the proposed TeV33 upgrade designed to yield an integrated 
luminosity of 25 fb~^. Our analysis is performed within the framework of the 
minimal supergravity model with gauge coupling unification and radiative 
electroweak symmetry breaking. For each of these three luminosity options, 
we delineate the regions of parameter space where jets plus missing energy plus 
0, 1, 2 (opposite sign and same-sign dileptons), and 3 isolated lepton signals 
from the cascade decays of sparticles should be visible above standard model 
backgrounds. We compare these with the parameter regions where signals in 
the clean isolated dilepton and trilepton channels (from chargino/neutralino 
and slepton production) should be observable. 



1 



Typeset using REVTgX 



I. INTRODUCTION 



The CDF and DO experiments have each accumulated an integrated luminosity of about 
100 pb~^ in Run I of the Fermilab Tevatron. An analysis of these data, which include 
elementary particle collisions at the highest energies accessible today, has already led to 
the discovery of the top quark and could reveal deviations from the expectations of 
the Standard Model (SM) of particle physics. These could take the form of new degrees of 
freedom (supersymmetry, technicolor, new gauge bosons, . . . ) or new effective interactions 
(quark compositeness, . . . ). At the very least, if no such deviation is observed, these data 
would serve to put phenomenological bounds on various extensions of the SM. Around 1999, 
the Main Injector (MI) is expected to begin operation: this should result in an order of 
magnitude increase in the size of the data sample. Tevatron experiments should then be 
sensitive to new physics processes with cross sections that are ten (three) times smaller 
than those accessible from an analysis of Run I data in rate-limited (background-limited) 
channels. 

Rather general arguments based on the instability of the SM's electroweak-breaking sec- 
tor to high mass scales suggest that it should break down at an energy scale A < 1 TeV, 
which is the raison d'etre for supercolliders such as the Large Hadron Collider (LHC) at 
CERN, or an e"'"e~ collider operating at a center of mass energy of ~1 TeV. If A is well 
below its upper bound, deviations from the SM may also manifest themselves at the Teva- 
tron, if sufficient integrated luminosity can be accumulated. Recently, the demise of the 
Superconducting Supercollider (SSC) has led several authors to propose that a luminos- 
ity upgrade beyond the Main Injector, designed to provide a data sample of 10 - 25 fb~^, 
could have a significantly improved reach for new physics well before the LHC commences 
operation. We will refer to this proposed upgrade as TeV33. We make no attempt here 
to assess the credibility of the TeV33 goals, either for the accelerator or for the existing or 
upgraded detectors. 

The purpose of this paper is to make a quantitative comparison between the capabilities 
of the current Tevatron, the MI upgrade, and the proposed TeV33 for the discovery of super- 
symmetry We focus on supersymmetry for several reasons. Supersymmetry (SUSY) 
provides the only known weakly coupled (and hence perturbatively calculable) framework 
that can naturally stabilize the Higgs sector of the SM. Since supersymmetry is a decou- 
pling theory, SUSY models reduce to the SM if sparticles are heavy: thus, SUSY models 
are at least as consistent as the SM when confronted with the precision data from LEP. 
SUSY models include a natural candidate for dark matter, and they can be consistently and 
simply embedded in a grand unified framework. A completely different reason for focusing 
on supersymmetry for the purpose of comparing the capabilities of different experimental 
facilities stems from the fact that SUSY models contain several new particles with a variety 
of quantum numbers: colored scalars and fermions, as well as corresponding colorless par- 
ticles. A supersymmetry skeptic could simply view our studies as providing a "theoretical 
laboratory" to compare the capabilities of various projected experimental facilities. 

Several studies of the SUSY reach of TeV33 already exist and results have been 

summarized in Ref. 0,11 . These studies do not all focus on the same SUSY reactions, differ 
significantly in the details of the computation of backgrounds (crucial in determining the 
reach), and present the final results in forms not amenable to direct comparison p|. In 
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view of the importance of this issue, we felt that a systematic study in which all signals are 
simultaneously studied using a common simulation would be useful in making an assessment 
of the increased capability of TeV33 over the already approved MI upgrade. Moreover, we 
present for the first time a comparison of the reach in various multilepton channels for the 
three luminosity options at the Tevatron. 

For definiteness, we work within the minimal supergravity framework, with assumptions 
(including grand unification and radiative electroweak symmetry breaking) detailed in our 
earlier studies [^|[ of SUSY signals at the Tevatron. The masses and couplings of all 
sparticles are then determined by just SUSY four parameters, 

• mo, the common scalar mass at the unification scale, 

• mi, the common gaugino mass at the unification scale, 

• Aq, the common soft SUSY breaking trilinear scalar coupling, and 

• tan f3, the ratio of the vacuum expectation values of the two Higgs fields, 

together with the sign of the Higgsino mass term /i. Our results also depend on mf, we 
shall take m^ = 170 GeV. The physical masses and couplings relevant to phenomenology 
are then obtained using the renormalization group equations. We emphasize that we 
use this framework for expediency. Indeed the assumptions underlying this framework may 
ultimately prove to be incorrect. This is largely irrelevant for our purpose, which is only to 
compare the SUSY reach of the different Tevatron upgrades. 

To orient the reader with the masses of various sparticles within this framework, we show, 
in Fig. 1, contours of the gluino (solid lines), squark (dashed lines) and lighter chargino 
(dotted lines) masses in the rriQ-mi plane. We fix = and illustrate the masses for 
a) tan /5 = 2, < 0, h) tan /3 = 2, > 0, c) tan /3 = 10, < 0, and d) tan /3 = 10, /i > 0. The 
gluino mass contours are not exactly horizontal because of the difference between the pole 



and running gluino masses |10]. The bricked regions are excluded by theoretical constraints 



detailed in Ref. [0] while the hatched region is excluded from the non-observation of any 



SUSY signal at the Tevatron [|TT|, or LEP [jT2[, and includes the recent mass limit |T3 



m,T, > 65 GeV from LEP1.5. 

The cascade decays of gluinos and squarks can lead to a variety of multijet plus multilep- 
ton event topologies. Contributions from associated production processes, while included, 
are known to be small. In addition, electroweak production of charginos and neutralinos can 
lead to hadronically quiet multilepton signals via which to search for supersymmetry. In this 
paper, we use ISA JET 7.16 |T^ to map out the supersymmetry reach in various channels at 
a 2 TeV pp collider, assuming an integrated luminosity of 

• 100 pb~^, roughly corresponding to what might be achieved from an analysis of the 
data from the current run; 

• 2 fb~^, corresponding to the approved MI upgrade; and 

• 25 fb~^, corresponding to what might be possible at the proposed TeV33 upgrade. 
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We assume identical detectors for each of these options. We neglect event pile-up effects, 
which have been shown to be small at least for the clean trilepton channel from W1Z2 
production 0. We also neglect any experimental difficulties associated with operating at 
the high luminosities implied by TeV33. Our purpose here is mainly to evaluate which 
regions of SUGRA parameter space can be explored via various discovery channels, and how 
these change depending on different Tevatron luminosity options. 

The remainder of this paper is organized as follows. In the next section, we briefly 
describe the event simulation that we use for the computation of various signals and back- 
grounds. In Sec. Ill, we focus on various multijet plus multilepton signals from the pro- 
duction and cascade decays of gluinos and squarks as well as from chargino and neutralino 
production, and map out the reach for supersymmetry in several of these channels for the 
three luminosity options introduced above. We present our results in the mo — mi plane for 
cases (a-d) introduced above. Variation with Aq is briefly addressed. Sec. IV is devoted to 
the study of the corresponding reach via electroweak production of charginos and neutrali- 
nos. We compare the total SUSY reach for the three Tevatron options and briefly discuss 
how these compare to the reach for other facilities such as the LHC or an e'^e~ linear collider 
operating at 500-1000 GeV in Sec. V. 



II. EVENT SIMULATION 

The implementation of the SUGRA framework into ISA JET has been described elsewhere 
and will not be repeated here. We generate all the lowest order 2 — 2 SUSY subprocesses 
in our simulation of n lepton plus m jet signals with m > 2 (except for unimportant s-channel 
Higgs boson mediated subprocesses). However, for the simulation of the clean multilepton 
signals, we have generated only slepton and chargino/neutralino events, since gluino and 
squark decays will very seldom yield final states without central jet activity. 

To model the experimental conditions at the Tevatron, we use the toy calorimeter sim- 
ulation package ISAPLT. We simulate calorimetry covering —4 < rj < 4 with cell size 
At) X A(j) = 0.1 X 0.0875. We take the hadronic (electromagnetic) energy resolution to be 
70%/ ^/E (15%/ VE). Jets are defined as hadronic clusters with Et > 15 GeV within a cone 
with AR = y/At]'^ + = 0.7. We require that \r]j\ < 3. Muons and electrons are classified 
as isolated if they have px > 7 GeV and \r){i)\ < 2.5 and if the visible activity within a 
cone of R = 0.3 about the lepton direction is less than £'5" (cone) = 5 GeV. In our analy- 
sis, we neglect multiple scattering effects and non-physics backgrounds from photon or jet 
misidentification, and we make no attempt 4to explicitly simulate any particular detector. 



III. REACH IN MULTILEPTON PLUS MULTIJET CHANNELS 

The multijet plus I^j, signal from gluino and squark production has been regarded as 
the most promising signature for supersymmetry at hadron colliders, provided that they 
are kinematically accessible. From a non-observation of any excess of events above SM 



expectations, the CDF and DO experiments [11 1 have already obtained the lower limits in 



the vicinity of 100-200 GeV on the masses of gluinos and squarks. These bounds have 
been derived within the framework of the minimal supersymmetric model, assuming ten 
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degenerate squark flavors. It has also been shown |T5|-p!^ that depending on rriq and rrig, 
other signatures with multilepton plus jets in the flnal state may also be observable with a 
data sample of C(IOO) pb~^ 

If gluinos and squarks are relatively heavy, electroweak production of charginos (l^i) 
and neutralinos {Z2) (which are expected to have masses ~ ^nig in models with a common 
gaugino mass at the uniflcation scale) may offer the best hope for SUSY detection at the 
Tevatron With the integrated luminosity that has been accumulated in Run I of the 
Tevatron, these signals should be on the verge of observability [ll9|-[2^. For an integrated 



luminosity in excess of 0(1 fb~ ) that should be available once the MI commences operation, 
the reach via the clean trilepton signal from the process pp —>■ W1Z2 — > iuZi + i'i'Zi may 
exceed that from gluino and squark reactions 1^-0 for a wide range of parameters. We will 
defer the discussion of these reactions to the next section, and focus our attention on the 
multijet plus multilepton signatures for now. 
We require, 

• jet multiplicity rijet > 2, 

• transverse sphericity St > 0.2, 

• > 40 GeV. 



As in our analysis of signals at the LHC p3[, we require an analysis cut, 

• EtUi), Et{]2) > E'j, and > E^, 

where the parameter is appropriately adjusted as described below. We further classify 
the events by their isolated lepton content as follows: 

• (A) I^rp events, with no isolated leptons. For this sample, we require that the missing 
energy not point along a jet, A0(^j,, Etj) > 30°; 

• (B) li events with exactly one isolated lepton with -E't(^) > 10 GeV. To reduce the 
background from W production, we also require Mr(£, > 100 GeV; 

• (C) Opposite sign (OS) dilepton events with exactly two unlike sign isolated leptons, 
where we require -Et(^i) > 10 GeV; 

• (D) Same sign (SS) dilepton events with exactly two same sign isolated leptons, again 
with Er(£i) > 10 GeV; 

• (E) 3i events, with exactly three isolated leptons with -Et(^i) > 10 GeV. We veto 
events with \M{i+i-) - Mz\ < 8 GeV. 

SUSY events are frequently rich in central b jets which can come from direct decays of 
gluinos and 6-squarks; from the decay Z2 bbZi, which can have an enhanced branching 



fraction pj]; or from the production of Higgs bosons in SUSY events. We have therefore 
studied the prospects for the detectability of tagged b signatures in the sample of and 1£ 
SUSY samples. In our analysis, we take the efficiency of tagging a 6-jet with Et > 15 GeV 
and \rih\ < 2 is 50%, and assume that the probability of mistagging other jets as a 6-jet is 
2%. We find that the signals for 6-tagged li events occur at unobservably small rates over 
most of the parameter space, so that for the most part, we will confine our attention to just 
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• (F) B + jets + I^rp events, where we veto any leptons, and require at least one tagged 
5-hadron. 

We have used ISAJET to compute these signals and the corresponding SM backgrounds 
which mainly come from, 

• W or Z production, in association with jets (additional leptons could arise from radi- 
ation), 

• ti production, 

• WW, WZ and ZZ pair production, where jets come from QCD radiation, and 

• QCD jet production, with I^j, due to mismeasurement of the jets. 

Some care must be exercised in the generation of the backgrounds. It is highly more 
likely that events with a hard initial scattering will pass the hard cuts that we have imposed 
to separate the SUSY signal from the background (especially for large values of Efj.). Since 
the cross section falls rapidly with the hard scattering p^, it is necessary to generate the 
backgrounds in various px bins to ensure that the relevant portions of the phase space are 
properly sampled. If this is not done, most of the events generated by the Monte Carlo pro- 
cedure will be too soft to pass the hard cuts, leading to an underestimate of the backgrounds. 



To avoid this we follow the procedure detailed in Ref . ||2^ , and generate background events 
in hard scattering pt in geometrically increasing bins between 25-400 GeV. 

The E^^ dependence of the various backgrounds is shown in Fig. 2 for (a) I^j, events, 
(b) li events, (c) OS dilepton events, (d) SS dilepton events, (e) 3i events, and (/) 6-|-jets-|- 
I^j. events. We use CTEQ2L structure functions throughout our analysis. The long- 
dashed line shows the background from ti production, while the QCD background is shown 
by the long-dashed-dotted line. The long-dashed double-dotted (triple-dotted) line shows 
the backgrounds from W + jets {Z -|- jets). Backgrounds from vector boson pair production 
are shown by the short-dashed (dotted) lines and are negligible, since the cross section for 
producing vector boson pairs together with at least two jets is rather small. The sum of all 
the backgrounds is shown as the solid line. We see from Fig. 2 that for all but the I^j, signal 
(and, to a lesser extent, the jets + 5 + -^^ signal), ti and 14^+jets backgrounds dominate over 
most of the range of For the signal, the background from Z ^ vv plus jet events is 



also significant and, in fact, dominates for very large values of E^. Backgrounds from 
QCD and vector boson pair production are negligible in this case. 

We have shown in our previous analyses of SUSY signals at the LHC that the reach 



can be optimized by adjusting the value of E'^ in the analyses. Events from relatively light 
gluinos and squarks tend to be softer but have larger cross sections than SUSY events when 
gluinos and squarks are quite heavy. While a modest value of E'^j, is probably optimal for 
the former case, the signal to background ratio as well as the statistical significance of the 
signal can be significantly improved by using larger values of E^ when attempting to extract 
a signal for very heavy sparticles. This should also be true of experiments at the Tevatron. 
As the size of the data sample increases, generally speaking, we may expect that we would 
obtain the maximum reach by increasing the value of E'j.. The analysis should, however, not 
be done for a single choice of this cut parameter that is optimized for the maximal reach 



6 



because this would result in a very low efficiency for signal detection if sparticles happened 



to be relatively light |23[. The optimal choice of E!j^ is also channel-dependent and somewhat 

sensitive to where we are in the parameter space. 

To underscore this, we have shown in Fig. 3 are total signal cross sections as a function of 
for the same event topologies {a)-{f) as in Fig. 2, but for several cases of input SUGRA 

parameters. In our illustration, we have fixed mi = 120 GeV so that rrig is roughly at 

350-380 GeV. In the ffist three cases, we choose tan/3 = 2, Aq = 0, fi <0, and take 

• (1) mo = 100 GeV (sohd) 

• (2) mo = 200 GeV (dashed) 

• (3) mo = 800 GeV (dotted) 

Increasing mo mainly increases the mean squark mass: in case 1, the decays g ^ qq are 
kinematically accessible, and of course, the production of squarks is substantial; in case 2, 
squarks of the first two generation are just heavier than the gluino, so that the dominant 
gluino decay is via g — >• bb; in case 3, all squarks are very heavy, and the gluino decays 
via three body modes. We have checked that flipping the sign of /i does not qualitatively 
change the results. In the remaining two cases 4 and 5, we fix mo = 200 GeV with /i > 0, 
and choose 

• (4) Ao = -400 GeV (dot-dashed) 

• (5) Ao = -430 GeV (dot-dot-dashed) 

These have been chosen so that the decay g tti is kinematically accessible: In case 4, 
ti decays via ti bWi, while in case 5 this channel is inaccessible and ti decays via the 



loop mode ti — > cZi p7| , |28[| . The size of the tagged 6-jet cross section is worth noting. 
The wiggles in the curve reflect the statistical fluctuations in the simulations. The total 
SM background is shown by the crosses. We see that if squarks are heavy (case 3, dotted) 
the background substantially exceeds the signal in all the channels and for the entire 
range where the signal cross section remains observable, and it appears unlikely that gluinos 
(whose mass is 380 GeV) will be detectable even at TeV33. We also see that for some signals 
{e.g. the OS channel), the choice of E!j^ significantly changes the signal to background ratio 
depending on the decay pattern of the gluino (compare cases 1 and 2 with 4 and 5). Clearly 
some optimization should be possible if the very large data sample as envisioned at TeV33 
become available. 

In order to assess the observability of the various SUSY signals at the Tevatron for 
integrated luminosities of 0.1 fb~^, 2 fb~^ and 25 fb~^, we have generated SUSY events 
for a grid of points in the mo — mi plane, for tan/3 = 2 and tan/3 = 10 for both signs 
of fi with Aq being fixed at zero. Our scan should be taken as spanning a representative 
portion of the parameter space for small and modest values of tan/3 [^. Since Aq mainly 



affects the phenomenology of the third generation, we expect the signals to be relatively 
insensitive to variation of Aq (except for signals involving third generation fermions, such 
as the cross section for events with tagged b jets. We will return to this point later. We 
have then run the generated SUSY events through our toy detector simulation and classified 
them into the various topologies A-F introduced above. Finally, for each SUGRA point, we 
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have checked whether the SUSY signals are observable above SM backgrounds for each of 
the three values of integrated luminosities. Here, we consider a signal to be observable if, 
for the given integrated luminosity, we have {i) at least 5 signal events, (m) the statistical 
significance of the signal exceeds "Scr", i.e. A^signai > 5i/^back, and (iii) iVgignai > 0.2A^back- 
The third requirement, while somewhat arbitrary, is to ensure that we do not classify the 
signal with = 25 GeV (see Fig. 3a) with a cross section of 300 fb to be observable above 
the background of 3x10"^ fb. It seems evident to us that because there is no characteristic 
kinematic distribution whose shape the signal would grossly distort, this "1% excess" , while 
"observable" according to criteria {i) and (ii), would be impossible to detect. We check the 
observability of SUSY events for = 15, 40, 60, 80, 100, 120 and 140 GeV and consider the 
signal to be observable if it is so for any one of the values of Ej^. 

The results of our computation are illustrated in Fig. 4-Fig. 9, for each of the event 
topologies A-F, respectively for the same four choices of other SUGRA parameters as 
in Fig. 1. The theoretically and experimentally excluded regions, denoted by bricks and 
hatches, are identical to Fig. 1. The regions of the mo — mi plane where a signal is ob- 
servable are shown for an integrated luminosity of 100 pb~^ (black squares) corresponding 
to the data sample of Run I, of 2 fb~^ (gray squares) corresponding to a data sample that 
might be available at the MI, and finally, of 25 fb^^ (white squares) corresponding to a data 
sample that might become available after a few years of operation of the proposed TeV33 
upgrade. 

Several comments are worth noting. 

• By comparing the distribution of the solid squares in Fig. 4-Fig. 9, we conclude that 
the channel should provide the maximal reach when the data of Run I is analyzed. 
Values of mi up to about 100 GeV should be accessible if mo < 200 GeV. For positive 

2 

values of // this region would already have been accessible via the chargino search 
at LEP1.5 (from Fig. 1 it is clear that Wi tends to be heavier for /i < 0, with all 
other parameters the same) which is why we have no solid squares in Fig. 4b and 
Fig. id. We stress, however, that this depends crucially on the assumed gaugino mass 
unification, and a direct for the gluino signal (even for parameters in the hatched 
regions) is extremely important. Also, note from Fig. 4a that Tevatron experiments 
may probe regions of parameter space not accessible at LEP2. 

• It is interesting to see that even with an integrated luminosity of 100 pb"-*^, we see that 
there are regions of the parameter space where there should be confirmatory signal 
also in other channels. Except for isolated corners in the small mo region (where the 
leptonic branching fraction of charginos and neutralinos may be enhanced because 
their two-body leptonic decays are accessible), the leptonic or tagged-6 channels are 
unlikely to be discovery channels for supersymmetry at least for Run I of the Tevatron. 

• With a data sample of 2 fb~^ Tevatron experiments should be able to probe mi 

2 

values up to 150 GeV (corresponding to rrig ~ 400 GeV) in the l^j, channel if mo ^ 
200 GeV. This is considerably beyond the reach of LEP2. At the high values of mi 
that should be explorable at the MI, it is worth noting that electroweak production 
of charginos and neutralinos is a significant contributor to even the channel: the 
production of gluinos and squarks is kinematically suppressed, so that the signals 
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from the clcctroweak production of charginos and neutralinos which have masses ~ 
(| — ^)mg become independent. The same is true in the leptonic channels — in the 
3£ channel, the jets then mainly come from QCD radiation, which is included in the 
shower approximation in ISA JET. For larger values of mo, squarks become inaccessible, 
and the range of mi that might be probed becomes smaller. Not only is the reach 
at the MI considerably larger, we see that there should be confirmatory signals in the 
various leptonic channels (particularly if mo is not very large) and even in the tagged 
6-channel. Once again, we see that the maximum reach is obtained in the channel, 
and that the trilepton channel enhances the reach only for isolated values of SUGRA 
parameters. 

• We see from Figs. 4 that TeV33 should be able to probe mi values about 25 GeV 
beyond what may be explored at the MI in the I^j, channel; this corresponds to an 
increase in the gluino mass reach of about 65 GeV. The increase in the reach may be 
somewhat larger in the leptonic channels, as can be seen from Figs. 5-8. However, the 
maximal reach is still obtained via the channel. We thus conclude that if Tcvatron 
experiments are able to accumulate about 25 fb~^ of data, they would be able to probe 
gluino masses about 20-25% beyond what might be probed at the MI {rrig ~ 300 — 500 
GeV). It is, however, important to note that with such a large data sample, one would 
be almost guaranteed to see a signal in several channels over much of the parameter 
range where there is a signal in the channel. 

• Turning to the prospects for identifying tagged 6-jets in the sample of events, we 
see from Fig. 9 that it is unlikely from an analysis of the data from Run I. 6-tagged I^j, 
events should be identifiable over a significant range of parameters at both the MI or 
TeV33. Measurement of the tagged b cross section, as well as the multiplicity of 6-jets 
in SUSY events, could serve to yield information about the A-parameter or the 6-quark 
Yukawa coupling, if a sufficient number of SUSY events can be accumulated. We warn 
the reader that the observability of this signal is rather sensitive to the assumptions 
about 6-tagging and mistagging. Wc have checked that the signal to background ratio 
becomes smaller if the mistagging probability is taken to be zero instead of 2%; i.e. 
the signal contains a larger fraction of mistagged events than the background. This is 
mainly due to the fact that the main SM background to the fe-tagged sample comes 
from top quark production, and the background rate from misidentification of QCD 
jets in W and Z events adds little to the main physics background. In addition, the 
larger jet multiplicity, and the resulting greater probability of mistagging jets in the 
signal sample should lead to a larger "fake" background for the signal events relative 
to SM events. 

Up to now, we have fixed the soft breaking parameter Aq — 0. In order to investigate 
the sensitivity of the cross sections to variation of Aq we show the Tevatron cross sections 
after cuts in the six channels A-F in Fig. 10. We fix mo = 100 GeV and mi = 120 GeV, 
which yields rriq ~ rrig — 30 GeV ^315 GeV for Aq = 0. The masses of the gluino and the 
first two generations of squarks are only very weakly dependent on Aq. Thus, for the cases 
studied in Fig. 10, qq, qg and gg processes all contribute to the signals, and gluinos decay 
to squarks. We examine values of for which the squared stop masses are positive. The 
following features are worth noting: 
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1. In case (a) the various topological cross sections are quite insensitive to Aq. As we 
vary Aq over the complete range, sparticle masses do not change very much, and no 
thresholds for new decays are encountered. The gluinos and squarks each decay via 
two body modes, with branching fractions relatively insensitive to Aq. 

2. In case (6), as we consider decreasing values of Aq, starting with large positive values, 
the cross sections are relatively constant until Aq becomes negative enough so that 
ti, the lighter of the t-squarks, becomes accessible via gluino decays. The scalar top 
decays via ti bWi, so that the number of 6-jets in gluino events is significantly 
increased. Direct production of titi pairs also contributes to the increase in the tagged 



6-jet cross section and should be independently detectable if m^^ is not too large. 
The increased cross sections in the OS and SS lepton channels also have their origin in 
g —>■ til decays followed by the semileptonic decays of the top family; iiii production 
contributes to the increase in the OS as well as the ^j. cross sections. 

Cases (c) and (d) show roughly similar features as case (6) — as Aq decreases from large 
positive values, the cross sections again show an increase because ti becomes relatively 
light. However, a new feature is now seen for large negative values of Aq in case (c) — 
the cross section shows a sharp increase, while all other cross sections drop sharply. 
This occurs when ti becomes lighter than the chargino, so that the loop decay mode 
ii — >■ cZi dominates ii decays, leading to a drop in the 6-jet multiplicity from gluino 
and stop decays. The leptonic cross sections decrease because the chargino now decays 
via Wi hti with a branching fraction of essentially 100%; because — m^-J is 

small, the 6-jets are presumably too soft to satisfy the tagging requirements. If SUSY 
parameters happen to be in this range, the promising clean trilepton signal from W1Z2 
production will be absent, and events will be the main signature for chargino pair 
production. Finally, we remark that in case [d), {m^^ — m^-J is somewhat larger (at 
extreme negative values of Aq) and the 6-jet appears to be sufficiently hard. As a 
result, while the leptonic cross sections exhibit a sharp decrease, the I^j. and tagged b 
topologies do not. 



IV. REACH VIA CLEAN MULTILEPTON CHANNELS 

If gluinos and squarks are beyond the reach of the Tevatron, it is possible that SUSY 
might manifest itself through electroweak processes, via events with two or more hard 
isolated leptons and I^j. but with essentially no jet activity. We will refer to these as 
clean multilepton channels. The spectacular trilepton event signature from the reaction 
pp W1Z2 — * iuZi + i'i'Zi has recently been the focus of considerable attention but 
observable signals may also be present in dilepton channels 0. Sources of clean dilepton 
events include chargino and slepton pair production. Since we have studied these channels 
in detail elsewhere, here we will mostly focus on the comparison of the capabilities of the 
three luminosity options, and refer the reader to our earlier paper ^ for details about the 
features of the SUSY signal and SM backgrounds. 
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A. Trilepton Channel 



To facilitate efficient computation of the signal levels in the clean 3i channel we generated 
all chargino, neutralino and slepton production subprocesses using ISAJET even though the 
signal is dominated by qq W1Z2 production. Gluino and squark production very seldom 
leads to events without any jet activity. Following our earlier analysis we implement the 
following cuts: 

• We require 3 isolated leptons within \rii\ < 2.5 in each event, with prih) > 20 GeV, 
Pt(^2) > 15 GeV, and pt(4) > 10 GeV; 

• We require I^j, > 25 GeV; 

• We require that the invariant mass of any opposite-sign, same flavor dilepton pair not 
reconstruct the Z mass, i.e. we require that \m{ii) — Mz\ > 10 GeV; 

• Finally require the events to be clean, i.e. we veto events with jets. 

SM backgrounds to the clean 3£ sample mainly come from ti and W Z production. Lepton 
isolation plays a crucial role in reducing the top quark background. After the cuts, we are 
left with [^] a SM physics background of just 0.2 fb (mainly from WZ events where Z — >■ rf 
with subsequent leptonic decays of the rs). Assuming that detector-dependent backgrounds 
from particle misidentification or jets or photons faking a lepton are under control, the 
observation of just a handful of such events in the current run (or at the MI) could signal 
new physics. 

In Fig. 11, we show the regions of the tuq — mi plane where the SUSY signal is expected 

2 

to be observable at the 5(T level (with a minimum of 5 signal events) at the MI (gray 
squares) and at TeV33 (white squares) for the same cases (a)-(c?) as in Fig. 1. We did 
not find any region (compatible with current experimental constraints) where this signal 
might be observable for an integrated luminosity of 100 pb~^ only because our cuts in this 
channel were optimized for the higher Tevatron luminosity options. The following features 
are worthy of note. 

1. We see that for small values of mo (for which sleptons are light enough so that two body 
slepton-lepton decays of charginos and neutralinos are kinematically allowed, while qq 
modes are forbidden) the MI reach extends to mi = 200 — 250 GeV depending on 

tan/3 and the sign of /i; the gap around mo — 50 GeV is where the decay Z2 IlC- 
becomes kinematically forbidden and where Z2 dominantly decays via Z2 vv- For 
still larger values of mo sneutrinos become too heavy and Z2 ini dominates, so that 
the trilepton signal is again observable. The kinematic boundaries for these various 
two body slepton decay channels of Z2 are shown by the three slanting contours in the 
small mo region of the Figure. Ultimately, of course, the sleptons become virtual, so 
that the neutralino decays via three-body modes. Then the hadronic decay of Z2 are 
no longer negligible, leading to a reduction in the trilepton cross section. In fact for 
mo > 150 — 200 GeV, the 3£ signal falls below the observable level at the MI except 
in the small tan/3, /i < case (a), where the range of mi that can be probed slowly 
decreases with increasing mo. 
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2. TeV33 should be able to probe the SUSY signal in this channel for substantial regions 
of parameter space not accessible at the MI. This should not be surprising, since the 
signal is essentially rate limited, assuming that non-physics instrumental backgrounds 
will be in control even in the high luminosity environment. Experiments at TeV33 
should not only be able to fill in most of the "gaps" at small mo where the signal is 
not observable at the MI, but should also be able to substantially extend the mo range 
where the signal might be observable. 

3. At TeV33 the clean trilepton signal may be observable well beyond where the "spoiler" 
decay modes of Z2 become accessible if two body slepton decays of Z2 are kinematically 
accessible; the kinematic boundaries for these decays are shown by the approximately 
horizontal contours. For a limited range of parameters, the TeV33 reach extends out 
to mi = 280 GeV which corresponds to vtig = 700 GeV! There are, however, even 
larger ranges of parameters where this signal will not be observable at TeV33 even if 
the chargino mass is at its current experimental bound. This has been traced to a 
negative interference term which leads to a large suppression of the leptonic decay of 

1^,0,^]. It is for the same reason that the signal is not observable all the way up 



to the limit of the "spoiler" mode in cases {b)-{d). The non-observation of a trilepton 
signal at TeV33 cannot, therefore, be translated into a lower limit on the Wi and Z2 
masses. 



B. Dilepton Channels 

We have also investigated the possibility of discovering supersymmetry in the OS clean 
acoUinear dilepton channel. Although the main contributions to the signal might be ex- 
pected to come from chargino and slepton pair production, we generate all possible SUSY 
production reactions and implement the following set of cuts, designed to extract signal from 
these backgrounds 0. 

• We require exactly two isolated OS (either e or /i ) leptons in each event, with pri^i) > 
10 GeV and prih) > 7 GeV, and \ri{£)\ < 2.5. In addition, we require no jets, which 
effectively reduces most of the ti background. 

• We require I^j, > 25 GeV to remove backgrounds from Drell-Yan dilepton production, 
and also the bulk of the background from 7* , Z ^ rf decay. 

• We require (f){ii) < 150°, to further reduce 7*, Z ^ rf background. 

• We require the Z mass cut: invariant mass of any opposite-sign, same flavor dilepton 
pair not reconstruct the Z mass, i.e. m{ii) 7^ Mz ± 10 GeV. 

The SM background, after these cuts, is about 44 fb, bulk of which comes from WW 
production. The regions of parameter space where the signal might be observable in this 
channel is shown in Fig. 12, again for the same cases {a)-{d). The horizontal (inclined) 
contours denote the kinematic boundaries for the decay(s) Wi — >• WZi (Wi Ilu, utj. A 
comparison with Fig. 11 shows that this region is a subset of the regions that might be 
explored via the clean 3£ channel, both at the MI as well as at TeV33. In particular, the 
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clean 2i channel does not lead to an observable signal in the regions of parameter space 
where the 3i signal is unobservable due to the strong suppression of the branching fraction 
for leptonic Z2 decays [^. Also, for the most part, the region of parameters that might be 
probed at TeV33 via the clean dilepton channel may be probed at the MI via the clean 3i 
channel. 

It should be possible to confirm a signal from W1Z2 production by searching for OS like- 
flavor dileptons plus jets plus events which arise when the chargino decays hadronically 
and the neutralino leptonically. The invariant mass of the lepton pair in these events should 
match up with the mass of the OS like-flavor dilepton pair in the sample of clean 

trilepton events. These events generally tend to be softer than events from the production of 
gluinos and squarks studied in Sec. III. We have examined where in SUGRA parameter space 
this jetty dilepton signal might be observable, with the same cuts as in Ref. 0. Because 
these events are detectable only over a subset of parameters where the clean 3i signal is 
observable, this channel does not lead to an improved reach. We will content ourselves with 
a qualitative discussion of this channel. 

The SM background which, after cuts, mainly comes from vector boson pair production 
and Drell-Yan rf production is about 25 fb. The OS lepton plus jets plus I^rp signal region is 
observable over a subset of the regions where the trilepton signal is observable. This region 
has roughly similar shape as the region in Fig. 11, but covers only just over half as many 
points. For small values of ttiq where the charginos predominantly decay via the two-body 
lepton-slepton mode, these events mainly come from slepton and chargino pair production 
with jets coming from QCD radiation. For mo > 200 GeV, the signal is observable for a 
rather limited range of parameters except in case (a) where mi values up to 100 (140) GeV 
may be probed at the MI (TeV33). 

It is worth mentioning that a measurement of a{3i)/a{i'^i~ +iets + I^j,) could potentially 
yield a measure of the hadronic branching fraction of Wi, provided that these events could 
be separated from other SUSY sources (gluinos and squarks) which lead to the same event 
topology. For very small mg, we have already mentioned that WiWi, with jets coming from 
QCD radiation, is the biggest source of these events. We have checked that for larger values 
of mo the rij > 2 plus OS dilepton events, with the cuts of Ref. indeed come mainly from 
W1Z2 production: for example, for (mo,mi) = (300 GeV, 120 GeV) in case (a) almost 90% 

of these events originate in the chain Wi —>■ qqZi,Z2 —>■ £^£^Zi. Even for the (100 GeV, 
80 GeV) case, this fraction exceeds 50%; several other signals should be expected in such a 
scenario. Although we have not attempted to explore this, it may be interesting to examine 



whether the determination of the chargino decay pattern at a high luminosity upgrade ||33 



such as TeV33 is indeed viable. (At the LHC, one may be swamped by squark and gluino 
production processes.) 



V. SUMMARY AND CONCLUSIONS 



Motivated by the recent proposal for upgrading the luminosity of the Fermilab Teva- 
tron by another order of magnitude beyond the Main Injector era, we have examined the 
impact that experiments at such a facility (should its construction prove technically and fls- 
cally feasible) would have on the search for supersymmetric particles. We have used ISAJET 
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to compute the reach of TeV33, under the optimistic assumption of an integrated luminosity 
of 25 fb^^ and compared it with the corresponding reach of the MI, for which we assume a 
data sample of 2 fb~^. In our analysis, we simply assume identical detector performance at 
the two facilities. We make no representation as to whether upgrades of the detector and 
the electronics that would be necessary for funtioning in this high luminosity environment 
can be achieved in a timely fashion. Thus, while our estimates for the reach of the MI are 
probably realistic, our conclusions regarding the reach of TeV33 might be on the optimistic 
side. We have also discussed what might be possible from an analysis of the current Tevatron 
data sample of ~ 100 pb~^. This enables us to compare how experiments at the MI will 
improve on what we can learn from the data sample of Run I of the Tevatron, and put in 
perspective the capabilities of any upgrades that might be possible in the future. 

For definiteness, we have adopted the SUGRA framework (with its assumptions about 
universal scalar and gaugino masses at an ultra-high scale) for our analysis. We have exam- 
ined the multijets plus I^j, plus = 1-3 lepton channels (Sec. Ill) as well as the hadronically 
quiet dilepton and trilepton channels (Sec. IV) and delineated regions of the mQ-mi plane 
where the various signals might be observable above SM backgrounds for the three values 
of integrated luminosity mentioned above. Our main result is summarized in Fig. 13, where 
we show the regions of SUGRA parameter space where at least one of the SUSY signals is 
observable with the criteria defined above. We show our results for the same cases {a)-{d) 
in Fig. 1. The black squares are the points that can be probed at the "5(j" level with an 
integrated luminosity of 100 pb~^, while the gray (white) squares are where the signal should 
be observable with a luminosity of 2 fb~^ (25 fb~^). Several features are worth emphasizing: 

• The analysis of the data from the current run will allow experiments at the Tevatron 
to probe only a little beyond the current experimental bounds. The most promising 
channel is the multijets plus I^j, channel from the production of gluinos and squarks; 
beyond mi = 100 GeV {rrig ~ 300 GeV), their production cross section is kinematically 
suppressed, while the signals from electroweak chargino/neutralino production are 
still rate limited. We stress, however, that Tevatron experiments are direct probes of 
gluinos and squarks, and it is important to look for their signals even below the hatched 
region, since it is entirely possible that the assumption of gaugino mass unification 
(which is crucial to translate the LEP chargino mass bound to a bound in this plane) 
may prove to be incorrect. 

• Experiments at the MI should probe a significant portion of the rriQ-mi plane, as 
can be seen from the distribution of gray squares in Fig. 13. Here, the important 
contributing channels are the multijets plus Ifj, channel and, especially in the small 
mo region, the clean 3£ channel. We see that the region of the tuq — mi plane that 
can be explored is sensitive to both tan/3 and sgn//: for favourable vafues of these 
parameters, the experiments may probe mi as large as 250 GeV if mg < 150 GeV 
via the clean trilepton channel; but for somewhat larger values of mg, we see from 
Figs. 13 {h)-{d) that there may be no SUSY signal at the MI even if charginos are at 
their current experimental bound. 

• Assuming an integrated luminosity of 25 fb^-*^, we see that experiments at TeV33 may 
be able to substantially expand the region where a SUSY signal might be observable. 
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Generally speaking, the reach is most enhanced in the rate-limited 3£ channel (where 
the gain is proportional to the integrated luminosity (L), unlike the background-limited 
multi-jet channels where the statistical significance of the signal improves only as VL). 
Nevertheless, it is important to note that there are significant ranges of parameters 
which cannot be efficiently probed at TeV33 even if charginos are just beyond the reach 
of LEP2 (recall that LEP2 will probe the region ttl^^ < 80 — 90 GeV well before the 
TeV33 commences operation). It will thus be difficult to obtain an unambiguous lower 
bound on sparticle masses if no SUSY signal is observed in experiments at TeV33. 

• We also see from Fig. 13 that experiments at the MI and TeV33 will substantially 
expand the SUSY reach of the Tevatron. An important virtue of TeV33, however, is 
that there is an observable SUSY signal in several channels over a substantial portion 
of parameter space where SUSY should be detectable. Furthermore, if a signal for 
supersymmetry is found at the MI, then the order of magnitude increase in collider 
luminosity at TeV33 should allow more detailed information about the underlying 
SUSY parameters to be extracted from the event sample. 

Our analysis was performed within the minimal framework where the conservation of im- 
parity is assumed. If instead i?-parity is violated by baryon number non-conserving operators 
so that the LSP decays hadronically, the as well as the multilepton signals are greatly 
diminished, and the reach may be significantly smaller p^] than that outlined in Fig. 13. 



At this point it is worth recalling that in order for supersymmetry to ameliorate the 
fine-tuning problem of the SM, sparticles cannot be arbitrarily heavy 0: qualitatively, one 
requires that sparticles are not much heavier than the weak scale. Several authors p5|j36[| 
have attempted to quantify this and obtained upper limits on sparticle masses. While these 
bounds are admittedly subjective, they could be regarded as providing rough benchmarks 



for future facilities; e.g. Anderson and Castano ||3^ have argued that the most favoured 
region from this point of view is where mi ^ 150-200 GeV, mo ^ 200-300 GeV. It is also 
interesting to note that the lightest neutralino would be an acceptable mixed dark matter 



candidate if SUGRA parameters happen to be in this range |3^. These arguments suggest 



that experiments at the MI and TeV33 would probe some of the most promising regions 
of SUGRA parameter space. We believe, however, that while fine-tuning and cosmological 
considerations are indeed suggestive, the upper bounds on sparticle masses that are obtained 
from these should be regarded as qualitative, and that a sufficient "safety margin" should 
be allowed for any experiment that is designed to decisively confirm or exclude weak scale 
supersymmetry. This is only possible at hadron supercolliders such as the LHC where it is 
possible to probe the rriQ-mi plane over the region mo < 1.5 TeV, mi < 800 GeV |2^, or at 



future electron-positron colliders operating at ^/s > 1-1.5 TeV, where it should be possible 
to probe chargino and slepton masses up to the beam energy. 
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FIGURES 



FIG. 1. Contours of squark (dashed), gluino (solid) and chargino (dotted) masses in the rriQ-mi 

2 

plane of the minimal SUGRA model. Frames are shown for a) tan/3 = 2, fi < 0, b) tan/? = 2, /U > 0, 
c) tan/? =10, /i < 0, and d) tan/? =10, // > 0. We take mt = 170 GeV and Aq = 0. The bricked 
regions are excluded by theoretical constraints discussed in Ref. while the shaded regions are 
excluded by experiment. 

FIG. 2. Component SM background cross sections as a function of the cut parameter 
defined in the text for multijet plus a) -^y, b) 1£ + I^j^, c) OS dileptons+^^i, d) SS dileptons+^jn, 
e) 3£ + I^rp and f) tagged B + I^j, event topologies with cuts as defined in Sec. Ill of the text. 
The backgrounds that we have computed are ti (long-dashed), QCD (long-dashed-single-dotted), 
W + jets (long-dashed-double-dotted), Z + jets (long-dashed-triple-dotted), WW (short-dashed), 
WZ (short-dashed-dotted) and ZZ (short-dashed-double-dotted) . The sum of these backgrounds 
is shown as the solid curve. 

FIG. 3. Variation of the SUSY signals and total SM background with Ej, for the same event 
toplogies as in Fig. 2. We have fixed mi = 120 GeV, tan /? = 2, and chosen Aq = Q and fi < with 
mo = 100 GeV (sohd), mo = 200 GeV (dashed) and mo = 800 GeV (dotted). In the other two 
cases, we take /j, > and fix ^o = —400 GeV (dot-dashed) and Aq = —430 GeV (dot-dot-dashed). 
The SM background level is denoted by crosses. 

FIG. 4. Regions of the mo — mi plane where the multijets plus Erp signal is observable at a 

_ 2 

2 TeV pp collider according to the criteria discussed in Sec. Ill of the text for the same choices 
of parameters as in Fig. 1. We consider three values for the integrated luminosity: 100 pb~^ 
corresponding to Run I of the Tevatron (black squares), 2 fb~^ which is expected to be accumulated 
at the Main Injector (gray squares) and 25 fb~^ to be accumulated at the proposed TeV33 (white 
squares). The bricked and hatched regions are the same as in Fig. 1. 

FIG. 5. The same as Fig. 4 except for the multijet plus li + channel. 
FIG. 6. The same as Fig. 4 except for the multijet plus OS dilepton +I^rp channel. 
FIG. 7. The same as Fig. 4 except for the multijet plus SS dilepton channel. 

FIG. 8. The same as Fig. 4 except for the multijet plus 3i + I^j, channel. 

FIG. 9. The same as Fig. 4 except for the multijet plus tagged B + I^j, channel. 

FIG. 10. The Aq dependence of the SUSY signal cross sections in the six channels in Fig. 3. 
for mo = 100 GeV, mi = 120 GeV, = 15 GeV for tan j3 = 2 and 10 and for either sign of /i. 
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FIG. 11. The same as Fig. 4 except for the clean, i.e. jet-free 3i + Ifrp channel discussed in 
Sec. IV. The three slanting contours (from left to right) mark the boundaries of the regions where 
the two body decays Z2 — > and IrI become forbidden, whereas the (almost) horizontal 

contours mark the corresponding boundaries for the "spoiler" decays Z2 ZZ\ or RiZ\ . Outside 
the hatched region, we found no points where this signal would be observable from an analysis of 
the data from Run I of the Tevatron. 

FIG. 12. The same as Fig. 11 except for the clean dilepton channel. The slanting contours, 
from left to right, mark the kinematic boundary for the decays W\ — > ii^v^vt, while the roughly 
horizontal line marks the boundary for the decay W\ WZ\. 

FIG. 13. The cumulative reach for supersymmetry of the Tevatron and its upgrades via any of 
the channels in Figs. 4-12 with the same labelling as in Fig. 4. 
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